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Abstract. We introduce a two-band Kondo-lattice model to describe 
ferromagnetic half-metals with local magnetic moments. In a model study, 
the electronic and magnetic properties are presented by temperature dependent 
magnetization curves, band-structures, spin polarizations and plasma frequencies. 
These are obtained from numerically evaluated equations, based on the single- 
electron Green functions. We show that the mutual influence between the 
itinerant electrons and the local magnetic moments is responsible for several 
phase transitions of the half-metals, namely first and second order magnetic phase 
transitions, as well as half-metal to semiconductor and half-metal to semimetal 
transitions. 
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1. Introduction 

Half- metals will play the key role in spintronics jlj, as do semiconductors in current 
electronics. Differently to semiconductors however, the Fermi level in ideal half-metals 
lies in a gap between valence and conduction bands only for one spin direction, whereas 
for the other spin direction the bands overlap. Thus the electrical current is exclusively 
due to spin polarized charge carriers, with promising consequences when exploiting 
the spin degree of freedom in logical devices. In semimetals, valence and conduction 
bands overlap regardless of the spin direction, leading to unpolarized currents. 

The investigation of half-metallic ferromagnets started with density functional 
theory calculations by de Groot et al. '5] on the Heusler alloy NiMnSb. Subsequently, 
other materials were identified as half-metals, as e.g. other Heusler alloys such as 
{Pt, Fe,Co)MnSb [3], ferromagnetic oxides as Cr02^ and Fe304 [5], and colossal 
magnetoresistance (CMR) systems like Lai-xSrxMnOs f6]. 

There are no pure elements which are half-metallic ferromagnets. The classical 
ferromagnets Co and Ni have fully polarized 3d-states at the Fermi edge, but there are 
also 4s-states at Ep preventing a fully polarized current. A mechanism must be found 
either to push the bottom of the 4s-Band above the Fermi edge or to press the latter 
below the band bottom. This is done by alloying or by forming an oxidic compound. 
Thus all known half-metallic ferromagnets contain more than one element. 

An important issue is the explanation of the gap, the origin of which is as equally 
diverse as the origin of half-metalicity. First-principles bandstructure calculations 
[H [7] predict such a gap at T = 0. What is the physical reason for it? What happens 
at finite temperature? Magnon and phonon effects may give rise to a depolarization 
so that strict half-metalicity appears to be limited to T = 0? 

In this paper we would like to model a local-moment half-metal to understand the 
basic physics of the gap-behavior. Candidates for such local-moment half-metals are 
the CMR-manganites, diluted ferromagnetic semiconductors such as Gai-xMn^As 
[H [9], and even concentrated ferromagnetic semiconductors as EuS [9]. We will 
describe these materials by a two-band Kondo-lattice model. This model (also known 
as sf model) describes the mutual influence of two well-defined subsystems, localized 
magnetic moments stemming from a partially filled electron shell (e.g. 4/), and 
itinerant electrons in a partially filled energy band. Both subsystems are coupled by 
an on-site interband exchange interaction. One well-known consequence of this mutual 
influence is the carrier-induced (anti)ferromagnetism in the local-moment subsystem 
(RKKY; Rudermann-Kittel-Kasuya-Yosida) , and another is the striking temperature- 
dependence of the bandstates resulting, e.g., in a red shift of the band bottom upon 
cooling below Tn [TOlfTT]. 

To describe local-moment half-metals we apply a two-band Kondo-lattice model 
(the principle idea is depicted in figure [T]). The itinerant electrons belong to a 
conduction and a valence band. For simplicity both are considered s-bands. They 
are exchange coupled to a system of localized magnetic moments. The exchange 
coupling is ferromagnetic (J'^ > 0) for conduction electrons and antiferromagnetic 
(J" < 0) for valence electrons. We demonstrate how the temperature-dependence of 
the bandstates can lead to a half-metallic groundstate with a transition to semimetallic 
or semiconducting behavior with increasing temperature. In a mean-field picture this is 
simply due to the fact that upon cooling, the spin-up conduction band is energetically 
shifted downwards and the spin-up valence band upwards, due to interband exchange 
between the itinerant charge carriers and the localized moments. This results in an 
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Figure 1. Prom left to right: free system (J = 0) witli tlie electronic parameters: 
bandwidths and of valence and conduction bands, respectively and 

intrinsic band gap F; on-site sf coupling J" < and J'^ > of the valence and 
conduction electrons with the local-moments Si\ half-metal with spin polarized 
electrons at the Fermi edge and mean magnetization (5^) of the local- moments 



increased concentration of spin-up electrons. The opposite happens to the spin-down 
electrons. We investigate this important phenomena in terms of the local-moment 
magnetization and the electron spin polarization, and try to explain them via the 
special behavior of the quasiparticle-bandstructure. 

In previous studies [121 O [TH [151 [K] i single-band Kondo-lattice models were 
used to describe half-metals. In this paper we investigate the two-band Kondo-lattice 
model, we had already proposed in [T7] to study the electronic properties of EuBg. 
Now, we extend the analytical methods by the modified-RKKY approach JJJ, which 
allows us to calculate the magnetic properties self-consistently instead of considering 
the net magnetization as a parameter following the Brillouin function. Furthermore, 
the focus of this paper lies on the exploration of characteristics of the introduced 
model and the prediction of possible properties of local-moment half-metals. This is 
done in a general model study rather than a specific case study. 

In the next section we introduce the model and outline the theory we used for 
solving the underlying many-body problem. The third section is then for a detailed 
discussion and interpretation of the results, followed in the last chapter by some 
conclusions. 



2. Theory 

The proposed two-band Kondo-lattice model for half-metals describes the electronic 
and magnetic systems which mutually influence each other. The free electronic system 
(left picture in figure[T]) consists of two s-bands with widths W" and W of valence and 
conduction bands, respectively. Between these bands we consider an intrinsic band 
gap r . As discussed above, when switching on the magnetic interaction J between the 
spins of the electrons in the bands and the local magnetic moments Si (figure[Tl middle 
picture), the bands will shift spin-dependently and form the half- metal (figure [l] right 
picture). 



Two-band ferromagnetic Kondo-lattice model for local-moment half-metals 



4 



2.1. Hamiltonian 

We describe the two-band Kondo-lattice model for local-moment half-metals by the 
following Hamiltonian 



v.c 



b 



(1) 



The index b stands for valence (v) and conduction band (c), i and j are different 
lattice sites and (T =1 or J, is the spin index of the electrons, a^^ and a^^ are the 
annihilation and creation operators, respectively. The chemical potential fj, will 
be chosen to result in equal hole and electron concentrations. 

The first term of the Hamiltonian ([T]) is the kinetic part. T^j represent the hopping 
integrals of the electrons and are connected to the Bloch energies Sk in reciprocal 
space by Fourier transformation Tj'j = J^k SfeS''^ "^^'^ . For the Bloch energies 
we use the simple cubic, tight-binding energy dispersion [15] 

Efc = 2b j — [cosik^a) -f cos(fcj^a) -f cos(fc^a)] ^ , (2) 

where k is the wavevector, a the lattice constant and zi, — 5i,c — Stv the band sign 
operator. The dispersion relations ^ can be recognized in the rigidly shifted 
spin-down valence band and spin-up conduction band in the top-most panels of 
figure m 

The second term of the Hamiltonian ([T]) describes the on-site magnetic sf coupling 
between the spins of the itinerant band electrons and the spin of the local- 
moments Si. The coupling strength is J — J'^ — —J", as already discussed 
above. 



2. 2. Single electron Green functions 

For the solution of the many body problem the Green functions method is used. From 
the Hamiltonian ([1]) we can formally derive the single electron Green functions 

GiAE) = z — 7—-, (3) 

where b = v,c stands again for the valence and conduction bands, respectively. The 
single electron Green functions ([3]) arc written down by use of the electron self-energy 
Tj^{E). Further details about the derivation of the self-energy for the Kondo-lattice 
model can be found in [191 |2Q]. The interpolating self-energy approach 19J is sufficient 
for the case of half metals, since we are only dealing with low itinerant carrier densities. 
For the low occupied conduction band, the self-energy from [19j is used 

with the abbreviation aa = S{1 -I- 5") — z„{S^){l -\- Za{S^)), the spin sign operator 
Za = <5o-T ~ ^ctX, the mean magnetization (S^), the free propagator Gq{E) = 
^{E ^ £fc + 1^)^^ and what would be the mean field result for the self-energy 

= -y^z^s^'). 
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For the valence band we also have a low concentration of itinerant carriers, 
the holes. The self-energy from [12] needs some adaptation using the electron-hole 
symmetry. This was done previously in fTT' and leads to the electron self-energy for 
the valence band 

i:jE) - M„ + J ) (5) 

With the self-energies (j4]) and ([5]) the Green functions ([3]) for valence and 
conduction bands can now be determined. 



2.3. Electronic properties 

The knowledge of the Green functions allows the calculation of the spectral density 

SiAE) = -^IniGuAE) (6) 

for any wavevector k. Evaluated for wavevectors along the standard symmetry points 
of the first Brillouin zone Equation ([6]) leads to the quasiparticle band structure 
presented in the next section. 

Furthermore, we will present the spin polarization of the conduction electrons 

n-s — n\ 

P=-^ ^, 7 

n-\ +ni 

with the spin dependent electron density 

n. = Y.J f-(E)SLiE)dE (8) 

k 

and the plasma frequency 



The effective masses m* are calculated via 



^_(d^£GUEL 



dE^ka 



(9) 



(10) 



1 

denotes 



at the r-point k = (000). In Equation (g]) f-{E) = exp (^^^ -|- 1 

the Fermi function and in Equation (fTO|) E^^ — Sk — fJ. + ReG^g.(i?£^) represent the 
quasiparticle resonance energies. The intrinsic effective mass ttt ~ (Gh^) I {Wa^) for 
W = 3.5eT^ and a = equals half the free electron mass mg. 



2.^. Magnetic properties 

For a consistent model study of the half-metals, we are also interested in the magnetic 
properties. They are derived in the modified- RKKY approach [11], where the 
Kondo lattice Hamiltonian is mapped to an effective Heisenberg Hamiltonian and the 
magnetic properties are calculated self-consistently, using the single electron Green 
functions. For the proposed two-band Kondo-lattice Hamiltonian (|T]) we will briefly 
recall the major steps in analogy to the derivation in [llj . 
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The mapping of Equation ([T]) to an effective Heisenberg Hamiltonian is achieved 
by averaging out the itinerant s-electron degrees of freedom, but retaining the operator 
character of the local moment spin operators, yielding 

H,f ^ ^ Hff = - ^ JifS. ■ S, (11) 

In the averaging procedure the expectation value {o^k+q a'^'ka''^ calculated in the 
electronic subspace. This is achieved by the restricted Green function G^'^' , , 
for which the equation of motion can be written down exactly. The iteration of the 
equation of motion was possible but would lead to higher orders of spin-operator 
products. Therefore, in the first iteration the single electron Green functions ([3]) 
are used for decoupling. This leads to the effective Heisenberg exchange integrals as 
functionals of the Kondo-lattice single electron Green functions 



v.c 



b 



\j'f j /-(i?) ilm iG%,{E)Y^G%,{E)\ AE. (12) 



The Green functions ([3]) are here used in their real space representation G\j^{E) 



N ^ G\^{E)&^'^^' and G\jQ{E) stands for the Green function of the free system 
(J = 0). 

With the effective exchange integrals (|12l) the effective Heisenberg Hamiltonian 
PT|) can now be used to calculate the magnetic properties in a standard procedure. We 
define a Green function of the local spin operators {{S^;e°'^i ^-f))' '^hose equation of 
motion is decoupled using the Tyablikov- approximation. With the Callen method [21] 
we eventually arrive at an equation for the mean magnetization 

il + S + ip)ip^''+^ + jS - ip )il + (^)25+i 

(l-f (^)2S+l _ (p2 

where (p — N^^'^^ exp( ) — 1 is the mean magnon number and E{q) — 

2(5*^) (Jq£o ~ Jq^) magnon energies. The exchange integrals are here used in the 
reciprocal space representation J^^ = N^^ J2ij Jfj^e'^'^'^'' ■ 

Since the single electron Green functions ([3|) are functionals of the mean 
magnetization (|13p and vice versa, we get a closed system of equations which will 
be solved self-consistently. The results from those calculations are presented for one 
parameter set in the following section. 



3. Results 



The proposed two band Kondo lattice model ([T]) implies a strong mutual influence 
between the electronic and magnetic systems of the half-metals. On the one hand, 
the influence of the local magnetic moments on the itinerant carriers - electrons in the 
conduction band and holes in the valence band - leads to a magnetization dependent 
band occupation. On the other hand, the magnetic exchange coupling via the itinerant 
carriers (modified-RKKY; mRKKY) is a functional of the band occupation itself, thus 
is the mean magnetization. This mutual influence of the two subsystems is reflected 
in the different phase transitions which will be presented in this section. 

For our qualitative description of ferromagnetic half-metals we will reduce the 
parameter space through setting the band widths of the valence and conduction bands 
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Figure 2. (color online) Temperature dependent magnetization curves for 
different intrinsic band gaps F, fixed coupling strengths J"^ = — J" = O.SeV, band 
widths = = 3.5eV and local spin quantum number S = ^. The solid 
lines show second order phase transitions and the dashed line a first order phase 
transition at Tq- The intermediate regime is represented by the dash-dotted line 
which shows a first order phase transition close to Tq and one of second order at 
Tq- In the numeric calculations Equation H13| l was solved self consistently, with 
the effective Heisenberg exchange integrals I I12I I and the electron Green functions 



to W" = W = 3.5eV, the local spin quantum number to S — ^ and the exchange 
coupling strengths to J° = — J" = 0.5eV. The variation of the intrinsic gap parameter 
r will shift the center of masses of valence and conduction bands resulting in different 
regimes of band occupations, leading to the different behavior of the half-metals. We 
will vary the intrinsic gap exemplarily in the full range of possible ferromagnetism. 

3.1. Magnetization 

The magnetization curves in figure [5] are various in their shapes due to the dynamic 
band occupation. Two different regimes can be distinguished, namely those with 
magnetic phase transitions of first order (dashed line) and second order (solid lines). 
Furthermore, one parameter set (dash-dotted line) with a first as well as second order 
phase transition is shown, representing the intermediate range. The critical exponents 
of the magnetization at the second order phase transitions at the Curie temperatures 
Tc are mean-field hke (/3 = 0.5 ± 0.05). 

In all band overlap exists at low temperatures and the itinerant carriers 

establish ferromagnetism. If, for increasing temperature, the carrier concentration 
becomes too low to maintain ferromagnetism via the RKKY-mechanism, the 
magnetization drops down suddenly (first order phase transition - dashed line). If the 
carrier concentration is sufficient up to the Curie temperature, the phase transition is 
continuous (second order - solid lines). In the intermediate range (dash-dotted line) the 
carrier concentration in the bands is not big enough to hold up the RKKY-mechanism 
above a temperature of 285K, which leads to a sudden drop in the magnetization (first 
order transition). Opposite to the case with even lower band occupation (dashed 
line), the magnetization doesn't drop down to zero, since the band occupation at 




Figure 3. Bandstructure for ferromagnetic saturation at zero temperature 
(top), (S^) = (middle) and in the paramagnetic state (bottom) |22| . The 
left-hand side shows the spin-down bandstructure, the right-hand side spin- 
up. The parameter are bandwidths = = S.SeV, intrinsic band gap 
r = JS = 1.75eV, sf-coupling strengths J'^ = —J" = 0.5eV and local spin 
quantum number 5 = ^. The Fermi level lies at = OeV. The spectral density is 
plotted as grey scale for wavevectors along the standard symmetry points of the 
first Brillouin zone for simple cubic lattices (Equation dSl*). 



temperatures between 285K and 300K is sufficient to hold up a reduced magnetization. 
This leads to a the second order phase transition at Te- 
la the next section we will have a closer look at how exactly the band structure 
changes under the variation of the net magnetization. 

3.2. Bandstructure 

The quasiparticle spectral density can be measured in (inverse) photoemission 
experiments. In figure [3] we plot the spin dependent bandstructure as grey scale 
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of the spectral density ^ at wavevectors along the standard symmetry points of 
the first Brillouin zone of simple cubic lattices. The degree of blackening is a 
measure of the height of the spectral density, the width of the dispersion refers to the 
quasiparticle damping. A sharp deep-black line refers to long-living quasiparticles. 
The representation fits the bare line-shape of a respective (inverse) photoemission 
experiment. The bandstructure is presented for three different mean magnetizations 
which were fixed during the calculations [22]. The intrinsic gap T = l.75eV in figure [3] 
is chosen for an easier discrimination between valence and conduction bands. For the 
intrinsic gap parameters used in figure [2] the bandstructures would look qualitatively 
the same as in figure [3l only the bands would be rigidly shifted towards each other by 
the difference in the intrinsic band gap. 

Before we come to a detailed discussion of the dynamic band occupation, 
responsible for the different phase transitions, we will explain the general properties 
that can be seen from the band structures of the two-band Kondo-lattice model. 

The upper right panel in figure [3] shows the two spin up bands. The conduction 
band is the undeformed tight-binding band. At zero temperature the crystal is 
ferromagnetically saturated (all local- moments are aligned parallel). Thus no spin 
exchange processes between the spin-up conduction electrons and the local-moments 
are possible, and the band is only rigidly shifted by —^JS. The same holds for the 
spin-down valence band (top left panel), due to the more than half-filling of the valence 
band. 

The faded regions in the bandstructure result from scattering processes of the 
electrons with the local- moments. Due to the scattering, the lifetime of those 
states is finite. In the sharp regions the spectral densities are delta functions and 
the quasiparticles have infinite lifetime. They cannot take part in spin exchange 
processes with the local-moments because no corresponding states with opposite spin 
are available. Instead one can think of virtual scattering for those states. At zero 
temperature the corresponding quasi particles are dressed by virtual clouds of magnons 
called magnetic polarons. This solution for the Kondo lattice model can be found 
analytically and is part of the interpolating self energy approach [T^] used for Equation 
Q and dni). 

With increasing temperature and decreasing magnetization (middle and lower 
panels), the bands almost split into two fiattened sub-bands. The small overlap of 
spin-up valence and conduction bands decreases with temperature and turns into a 
considerable bandgap in the paramagnetic state (lower panels). 

As already mentioned above for different intrinsic band gaps, like those in figure[21 
the bands in figure[3]would be rigidly shift towards each other resulting in greater band 
overlaps. The band overlap still changes with magnetization but varies in a different 
range. The different ranges are responsible for the different magnetization curves seen 
in figure [H since the effective direct exchange coupling between the local-moments is 
mediated by the itinerant carriers (mRKKY- mechanism). 

For intrinsic band gaps less than O.leV, the band occupation is sufficient 
to maintain the mRKKY-mechanism in the whole temperature range. Thus the 
magnetization decreases continuously with increasing temperature, and the magnetic 
phase transition at the Curie-temperature is of second order (solid lines in figure [2|) . 

For intrinsic band gaps greater than O.leV the band overlap, existent for large 
mean magnetizations vanishes while the magnetization decreases. Thus the band 
occupation goes to zero and the mRKKY-mechanism cannot be upheld any longer. 
When reaching the critical temperature, the magnetization drops down suddenly and 
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we see a first order magnetic phase transition (dashed hnes in figure [2]). 

For intrinsic band gaps around 0.7eT^, the magnetization curve in figure shows 
a first as well as second order phase transition. In this case the band overlap is 
sufficient to hold up a mean magnetization in the whole temperature range but it 
is not continuous. At the first critical temperature the magnetization drops down 
rapidly to continue at a smaller value continously up to the Curie-temperature. 

3.3. Spin-polarization and electron density 

The important property of half- metals is the spin polarization. In figure[3^ we present 
the spin polarization P, calculated from Equation ([7]) and in figure HJa the electron 
density n = n| -I- n^, from Equation ([8]). The magnetization was calculated self- 
consistently as before. 

One can see in the case T — 1.05ey (first order magnetic phase transition), the 
spin polarization and electron density drop down suddenly at Tc- Here, the sample 
undergoes a transition from half-metal to semiconductor. 

For the cases F < 0.7eV with second order magnetic phase transitions we have 
a smoothly decreasing spin polarization. These samples undergo a transition from 
half-metals to semimetals. 

The case F = O.leV, where a first and second order magnetic phase transition is 
present, shows a corresponding drop of electron density and spin-polarization at the 
same critical temperatures. 

3.4. Plasma frequency 

The plasma frequency can be obtained from optical absorption/reflectivity 
measurements. Figure Ht shows the plasma frequencies ujp, calculated from Equation 
^ for the half-metals in discussion. It depends on the electron density as well as 
the effective mass, which both change with temperature in the introduced two-band 
Kondo-lattice model. Although, the electron density varies more, compared to the 
effective mass. Thus, the variation of the plasma frequency mainly depends on the 
temperature dependent electron density. Once more one can see the different behavior 
of the first and second order phase transitions as described above. 

4. Conclusion 

We introduced a two band Kondo lattice model which describes valence and conduction 
electrons that are on-site coupled to local magnetic moments. By propagating 
through the lattice, the electrons are responsible for an effective direct exchange 
interaction (modified-RKKY) between the local magnetic moments, which enables 
ferromagnetism. 

The mutual influence between the spin dependent band overlap of valence and 
conduction bands and the mean magnetization lead to different phase transitions. 
On one hand it is possible to describe first as well as second order magnetic phase 
transitions. On the other hand one can explain the spin-electronic phase transitions: 
half- metal to semimetal and half- metal to semiconductor. 

For future spintronic device applications one could think of two possible scenarios: 

• We saw that changing the intrinsic bandgap would change the Curie-temperatures 
and kinds of phase transitions. Applying e.g. hydrostatic pressure on a sample. 
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Figure 4. (color online) (a)temperature dependent spin-polarization P, (b) 
temperature dependent electron density n and (c) temperature dependent plasma 
frequency Wp for different intrinsic band gaps F, fixed coupling strengths J'^ = 
— J" = O.beV , band widths = W" = S.SeV and local spin quantum number 
5 = ^ . The dashed line shows a transition from half- metal to semiconductor and 
the other lines from half-motal to semimetal. 



thus changing the lattice constant in a smaU range and hence the intrinsic 
bandgap, could result in a variation of the electronic and magnetic properties. 
This would be especially interesting around the Curie-temperature. 

• A stronger impact on the sample would have the application of an external 
magnetic field at temperatures around the Curie-temperature. Depending on the 
specific properties of the sample, one can think of magnetic field-driven changes 
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from the semiconducting (semimetallic) state to a half-metallic state. 

In conclusion, first it is necessary to find appropriate materials for which the 
theory is suitable. The search for that should focus on compounds with the typical 
(ferromagnetic > 0, antiferromagnetic J" < 0) Kondo lattice rare earths, such as 
Ce, Eu. Although the theory was developed for periodic crystals with local-moments 
on every lattice site, it is surely possible to expand the model to randomly distributed 
local-moment compounds, which should then show similar features. This possible 
future work would be interesting for the diluted magnetic semiconductors (DMS). 
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